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ABSTRACT 

We report on a 95 ks Chandra observation of the TeV emitting High Mass X- 
ray Binary LSI+61°303, using the ACIS-S camera in Continuos Clocking mode to 
search for a possible X-ray pulsar in this system. The observation was performed while 
the compact object was passing from phase 0.94 to 0.98 in its orbit around the Be 
companion star (hence close to the apastron passage). We did not find any periodic or 
quasi-periodic signal (at this orbital phase) in a frequency range of 0.005 — 175 Hz. We 
derived an average pulsed fraction 3er upper limit for the presence of a periodic signal 
of ^ 10% (although this limit is strongly dependent on the frequency and the energy 
band), the deepest limit ever reached for this object. Furthermore, the source appears 
highly variable in flux and spectrum even in this very small orbital phase range, in 
particular we detect two flares, lasting thousands of seconds, with a very hard X-ray 
spectrum with respect to the average source spectral distribution. The X-ray pulsed 
fraction limits we derived are lower than the pulsed fraction of any isolated rotational- 
powered pulsar, in particular having a TeV counterpart. In this scenario most of the 
X-ray emission of LSI-|-61 o 303 should necessarily come from the interwind or inner- 
pulsar wind zone shock rather than from the magnetosphere of the putative pulsar. On 
the other hand, very low X-ray pulsed fractions are not unseen in accreting neutron 
star systems, although we cannot at all exclude the black hole nature of the hosted 
compact object, a pulsar with a beam pointing away from our line of sight or spinning 
faster than ~ 5.6 ms, nor that pulsations might have a transient appearance in only a 
small fraction of the orbit. Furthermore, we did not find evidence for the previously 
suggested extended X-ray emission. 
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1 INTRODUCTION 

High Mass X-ray binaries (HMXBs) are relatively young 
(< 10 7 year) systems composed of a massive OB- type star 
and a compact object, either a neutron star or a black 
hole. These systems are generally bright X-rays emitters, 
due to matter from the OB star accreting onto the com- 
pact object. The majority of the HMXBs are Be/X-ray bi- 
naries, and some tens of them are O main-sequence or super- 
giant X-ray binary systems. Accretion takes place through 
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either the strong stellar wind of the optical companion or 
Roche-lobe overflow, where matter flows via the inner La- 
grangian point and an accretion disk to the compact object 
(see Frank, King & Raine 1992). A powerful X-ray source 
(Lx ~ 10 35 — 10 38 ergs _1 ) is produced, which ionizes al- 
most the whole stellar wind. The X-ray continuum spectra 
of HMXBs are often described by a power law with photon 
index a ~ 1 — 2 (modified at higher energies by an expo- 
nential cutoff). A spectrum of this form can be produced by 
inverse Compton scattering of soft X-rays by hot electrons in 
the accretion column near the compact object, and a part of 
this emission is scattered by the stellar wind of the massive 
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Figure 1. LSI+61°303's geometry considering the orbital solu- 
tion of Casares et al. (2005). The phases for Inferior conjunction 
(INFC), Superior conjunction (SUPC), periastron (P), and apas- 
tron (A) are marked, and the inclination is not taken into account 
in the plot. The orbit (black solid line), the massive star (MS; in 
orange) and the equatorial Be disk (in yellow) are roughly to 
scale, assuming a disk radius 7 times larger than the stellar ra- 
dius. The blue thick line indicates the orbital phases spanned by 
our Chandra observation (adapted from Sierpowska-Bartosik & 
Torres 2010). 



companion. This might sometimes result in a further non- 
thermal spectral component, but with a different absorption 
column depending on the orbital phase of the system. Fur- 
thermore, in some HMXBs a soft excess at ~ 0.1 — 2 keV is 
detected, very common in systems hosting pulsars (Hickox 
et al 2004). 

On top of the continuum model, spectral lines are often 
present in these systems, neutral and ionized, such as Fe, Si, 
Mg, Ar, N, Ca, mainly produced in the stellar wind or in 
the accretion disk (if any) illuminated by the strong X-ray 
emission of the compact object (see e.g. Cottam et al. 2001; 
van der Meer et al. 2005). 

To date there are only four or five binary systems known 
to be emitting up to TeV energies, all of them of the HMXB 
class: PSRB1259-63 (Aharonian et al. 2005a), LSI +61° 303 
(Albert et al. 2006), LS 5039 (Aharonian et al. 2005b, 2006), 
CygX-1 (Albert et al. 2007) and possibly the newly sug- 
gested candidate HESS J0632+057 (Aharonian et al. 2007, 
Hinton et al. 2009). These sources are rather peculiar with 
respect to the HMXB class. The first identified TeV bi- 
nary system was PSRB1259-63, a 3.4 year period binary 
hosting a 48 ms radio pulsar in an eccentric orbit around 
a Be-star (Johnston et al. 1992). The TeV emission from 
this object is thought to be associated with the radio pulsar 
wind and its interaction with the radiation field and mate- 
rial around the Be-star, and shows variable un-pulsed X-ray 
(Chernyakova et al. 2006, 2009) and radio emission (John- 
ston et al. 1999, 2005), in particular close to the periastron 
passage due to the interaction between the pulsar and the 
Be disk. PSRB1259-63's X-ray emission is well described 
by an absorbed power-law (Nh — 5 x 10 21 cm -2 ; F ~1.2- 



1.8), with an unabsorbed 1-10 keV flux varying between 0.1- 
3xlO" 11 ergs" 1 cm" 2 (Chernyakova et al. 2006, 2009). 

LSI+61°303 and LS 5039 are both much closer bina- 
ries, with orbital periods of 26.5 and 3.9 days, respectively, 
hosting a very massive star (B and O types) and a compact 
object, the nature of which is still unknown for both sources 
(see Fig.[T]for a schematic view of LS I +61°303's orbit). Fur- 
thermore, LS I +61°303 and LS 5039 share the quality of i) 
being the only two known 7-ray binaries detected in the GeV 
band (Hartman et al. 1999; Abdo et al. 2009a, 2009b), ii) 
they both show periodic TeV emission modulated by their 
orbital motion (Aharonian et al. 2006; Albert et al. 2008, 
2009), iii) variable X-ray emission (Paredes et al. 2007; Sidoli 
et al. 2006; Bosch-Ramon et al. 2005, 2007) correlated with 
the orbital motion, and iv) they have a variable radio coun- 
terpart (see also below). Their X-ray emission is highly vari- 
able on a ks timescale, and is for both objects characterized 
by an absorbed power-law spectrum (LS I +61° 303: Nh = 
6 x 10 21 cm~ 2 and T ~1.4-1.9; LS 5039 : N H = 7 x 10 21 cm~ 2 
and F ~1.4-1.6), with a 1-10 keV unabsorbed flux varying 
between 0.4 - 1.5 and 0.5 - 1.2 x lO^ergs'W" 2 , for 
LSI+61°303 and LS5039, respectively (Sidoli et al. 2006; 
Esposito et al. 2008; Kishishita et al. 2009). 

As the nature of the compact object is unknown in these 
two systems, it is not clear if their emission (especially in the 
TeV band) is due to a relativistic outflow from a rotational 
powered pulsar (as the case of PSRB1259-63) or accretion 
onto a black hole or a neutron star which might e.g. ac- 
celerate particles through a relativistic jet. Recently, a Soft 
Gamma Repeater (SGR)-like burst detected from the di- 
rection of LS I +61°303 raised the possibility of this system 
hosting a magnetar, although this is so far a very contro- 
versial hypothesis (see de Pasquale et al. 2008; Barthelmy 
et al 2008; Dubus & Giebels 2008; Rea & Torres 2008, Ray 
et al. 2008; Munoz-Arjonilla et al. 2008). In the accretion 
scenario context, the recent evidence for TeV emission from 
the black hole binary CygX-1, is very exciting (Albert et 
al. 2007). However, the very high energy phenomenology of 
CygX-1 appears different from that of these other sources: 
it has been in fact detected in TeV just once, and during a 
flaring state. The three other sources, instead, present per- 
sistent TeV emission modulated by the orbital motion, and 
which dominates their radiative outputs. For the new can- 
didate TeV binary HESS J0632+057 , there is still only very 
little known (e.g. no orbital period detected yet) to make 
any solid comparison (Aharonian et al. 2007, Hinton et al. 
2009; Maier et al. 2009). 

For LS I +61 o 303, extended radio emitting structures at 
angular extensions of 0.01-0.05" have been reported (Massi 
et al. 2001, 2004). This discovery has earlier supported its 
microquasar interpretation, as a claim for a jet-like feature. 
However, a subsequent large VLBI campaign (Dhawan et 
al. 2006), discovered rapid changes during the orbit, in the 
orientation of what seems to be a cometary tail consistent 
with being the result of a pulsar wind more than a jet -like 
structure. In fact, if due to a jet, the changing morphology 
of the radio emission along the orbit would require a highly 
unstable jet, and the shape of this structure is not expected 
to be reproduced orbit after orbit as current results indicate 
(Albert et al. 2008). Similarly, the discovery of a jet-like ra- 
dio structure in LS 5039 , initially prompted a microquasar 
interpretation (Paredes et al. 2000). However, Dubus (2006) 
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Figure 2. Chandra background-subtracted lightcurve of LSI+61°303 binned at 100s, and in the 0.3-10 keV energy range. The orbital 
phase <f> is also reported. 



showed that this extended radio emission could also be inter- 
preted as the result of a pulsar wind, with no need to invoke 
a radio jet. Furthermore, recent VLBA observations showed 
that a microquasar scenario cannot easily explain the ob- 
served morphology changes of this radio extended feature 
(Ribo et al. 2008). 

At this point, one of the long-standing open questions 
of high energy astrophysics is: are these two peculiar binary 
systems accelerating charged particles until TeV energies, 
e.g. through the powerful jets of an accreting neutron star 
or black hole, as it might be the case of the TeV flares from 
Cyg X-l? Or is the TeV emission of these binaries instead re- 
lated to a young and energetic rotational powered-pulsar as 
the case of PSRB1259-63? To try to address this question 
we performed deep Chandra observations of LSI+61°303 
(this paper) and LS 5039 (Rea et al. 2010, in preparation) 
in search for pulsed X-ray emission with very low pulsed 
fraction (e.g. because contaminated by the wind-shock emis- 
sion), which due to the higher background of previous X-ray 
observations might have been missed. In ijUwe report on the 
LS I +61°303 observation and analysis, and we present our 
results in followed by a discussion O- 



2 OBSERVATION AND DATA ANALYSIS 

The Advanced CCD Imaging Spectrometer (ACIS) cam- 
era on board of Chandra observed LSI+61°303 on 2008 
November 14th (start time 11:09:58 UT; Obs-ID 10052) for 
an exposure time of 95.668 ks in Continuos Clocking (CC) 
mode (FAINT). The CC observing mode provides a time 
resolution of 2.85 ms and imaging along a single direction. 
The source was positioned in the back-illuminated ACIS-S3 
CCD at the nominal target position. Standard processing 
of the data was performed by the Chandra X-ray Center to 
level 1 and level 2 (processing software DS ver. 7.6.11.9). 
The data were reprocessed using the CIAO software (ver. 
4.1) and the Chandra calibration files (CALDB ver. 4.2.0). 

Since in the CC mode the events are tagged with the 
times of the frame store, we corrected the times for the vari- 
able delay due to the spacecraft dithering and telescope flex- 
ure, starting from level 1 data and assuming that all photons 
were originally detected at the target position. Furthermore, 
data were filtered to exclude hot pixels, bad columns, and 




Figure 3. LS I +61°303 lightcurve in the 0.3-2 keV and 4-10 keV 
energy bands (top and middle rows, respectively), and their ratio 
(bottom row), all reported with a 1 ks bin time. Vertical lines 
separate the 5 intervals used in the analysis (of which we report 
the photon index (r) values). 



possible afterglow events (residual charge from the interac- 
tion of a cosmic ray in the CCD). Photon arrival times are 
in TDB and were referred to the barycenter of the Solar 
System using the JPL-DE405 ephemeris. 

In order to carry out a timing analysis, we extracted 
the events in the 0.3-10 keV energy range from a region 
of 5x5 pixels around the source position (RA 02:40:31.670, 
Dec +61:13:45.11) to reduce the background contamination, 
while the source spectrum was extracted from a rectangu- 
lar region of 5x25 pixels around the source position, and 
the background was taken independently from a source-free 
region in the same chip. We extracted the response matrix 
files (RMFs) and ancillary response files (ARFs), first cre- 
ating a weighted image, re-binning by a factor of 8, then 
using it to build the RMF file using the mkacisrmf tool, 
with an energy grid ranging from 0.3 to 10 keV in 5 eV in- 
crements. Using this RMF and the aspect histogram created 
with the aspect solution for this observation (asphist), we 
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generated the appropriate ARF file for the source position. 
The source ACIS-S count rate in the 0.3-f keV energy band 
was 0.648±0.003 counts s _1 (the lightcurve is shown in Fig. [2] 
and Fig.©. 



3 RESULTS 

Using the orbital ephemeris of Casares et al. (2005), namely 
an orbital period of 26.4960±0.0028 days and the zero- 
phase assumed at To=2,443,366.775 JD (defined as the radio 
phase; Gregory 2002), we calculated that the Chandra ob- 
servation was performed while the neutron star was passing 
from phase 0.94 (at the start of the observation) to phase 
0.98 (when it ends). Even if considering other orbital solu- 
tions (Grundstrom et al. 2006; Aragona et al. 2009), this 
corresponds for all assumed solutions to a part of the orbit 
soon after the apastron passage (see Fig.[T]). 

3.1 Timing analysis: pulsation search and flares 

We searched for periodic and quasi periodic signals with un- 
precedented sensitivity, thanks to the extremely low ACIS-S 
background (see for previous limits: Sidoli et al. 2006; Har- 
rison et al. 20023; Torres, Sierpowska-Bartosik & Rea 2010). 
Given the length of our observation (95 ks) and the timing 
resolution of the CC-mode (0.00285 s), we could in principle 
search for periodic signals in the 10 - 175 Hz frequency 
range (namely between the Nyquist frequency and half of 
the frequency resolution of our Chandra data; van der Klis 
1989), however, our data were sensitive only in the 0.005 
- 175 Hz frequency range, which we will then use all over 
this study. We studied the source power spectra performing 
Fast Fourier Transforms (FFTs; see Fig.[4| using the Xronos 
analysis software, and we did not find any periodic or quasi- 
periodic signal neither considering the total exposure time 
nor in each of the five time intervals shown in Fig. [3] Fur- 
thermore, we performed the same search dividing the data 
in three energy bands (0.3-2, 2-4, and 4-10 keV), again not 
finding any significant signal. In particular, for com- 
puting reasons, we performed an average over 9 FFTs with 
a bin time of 0.00285 s (Fig. [4] top-left panel), resulting in 
4,194,304 bins (~ 12 ks) for each of the 9 averaged power 
spectra (the averaged power spectrum had a xt distribution 
with 18 d.o.f.). For the power spectra produced for the three 
different energy bands, given the lower number of counts, we 
could use a single FFT (their power spectrum had then a xl 
distribution with 2 d.o.f.) . We took into account the number 
of bins searched, and the different d.o.f. of the noise power 
distribution (Vaughan et al. 1994; Israel & Stella 1996) in 
calculating the 3u detection upper limits reported in Fig.|4] 
. At very low frequencies, most probably due to the over- 
all variability and flares present in the X-ray lightcurve (see 
Fig.©, there is evidence for the onset of a red- noise com- 
ponent. We accounted for the red-noise component in the 

1 Note that the iJXTE pulsed fraction limits reported by Harri- 
son et al. 2002 (~ 6%) are calculated from the total count rate 
without correcting for the cosmic and instrumental background, 
which if corrected increases substantially the upper limit on the 
detectable pulsed fraction to ~ 32%. 



calculation of the detection significance by using a smooth- 
ing window technique (see Israel & Stella 1996 for further 
details). 

We computed the 3<r upper limits on the sinusoid semi- 
amplitude pulsed fraction {PF), computed according to 
Vaughan et al. (1994) and Israel & Stella (1996), which 
ranges in the 0.3-10 keV energy band between PF <7-15% 
(see Fig.[4] bottom-left panel). In Fig. [4] (right column) we 
calculate the same pulsed fraction limits as a function of the 
energy band. It is straightforward that the lower number of 
counts present in the energy-selected lightcurves causes the 
energy-dependent PF limits to be slightly larger than those 
derived using the whole dataset. Note that given the long 
orbital period of this binary system (~ 26.5 days; Casares 
et al. 2005), we could reliably ignore the effect of such orbit 
during our pulsar searches, since our exposure time is much 
less than 10% of orbital period (Vaughan et al. 1994). 

We looked for the presence of SGR-like bursts by bin- 
ning the counts in intervals of 0.1s and searching for ex- 
cesses above a count threshold corresponding to a chance 
occurrence of 0.1% (taking into account the total number of 
bins), but we did not find any significant short burst. On 
the other hand, on top of a slow variability of a factor of 
~ 2, we significantly detect two flares lasting ~3 and 1.5 ks 
(not unseen in this source: Sidoli et al. 2006; Paredes et al. 
2007; Esposito et al. 2007; Smith et al. 2009) during the 
minimum of the lightcurve, with a rise-time of ~ 100 s, and 
peaking after about 28 and 33 ks from the observation start 
(see Fig.© and [3]). 



3.2 Spatial analysis: search for extended emission 

We studied the radial profile of LSI+61°303, using the CC 
mode one-dimensional imaging capability, to search for any 
evidence for extended X-ray emission (see also Paredes et al. 
2007). We first generated an image of the one-dimensional 
strip in the 0.3-10 keV band, and then subtracted the back- 
ground count rate to remove instrumental and cosmic X- 
ray background. We then produced a one-dimensional sur- 
face brightness distribution using the same method as for 
a two-dimensional radial profile, namely extracting photons 
from 125 annuli centered on the source position and of 2 
pixel wide (~0.99"). We then simulated a Chart/MARX 
Point Spread Function (PSF) and extracted photons from 
the same annular regions. We did not find evidence for a sig- 
nificant extended emission or structures within 30" around 
our target (see Fig.©. As a further check, given the possi- 
ble uncertainties of the one-dimensional CC-mode PSF, we 
extracted with the same method a similar profile from a 
source observed in the same mode, and known to be point- 
like, namely the AMCVn RX J0806+1527 (Obs-ID 4523; Is- 
rael et al. 2003). This source has a spectrum similar to our 
source, this is why we chose it as a good candidate for the 
PSF comparison. As shown in the inset of Fig. [5] we do not 
find any excess in the ratio between the surface brightness 
of LSI+61°303 and RX J0806+1527, but only a constant 
flux difference of ~ 31, which corresponds to the different 
brightness of these sources. 
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Figure 4. Left panels: Top row: power spectra of LSI+61°303 in the 0.3— lOkeV energy band binned at 0.00285 s with superimposed 
the 3<t upper limits on the detection of a periodic signal on top of the Poisson noise (we used here logarithmic scales to better show the 
low frequency range). Bottom row: pulsed fraction 3<r upper limits on the semi-amplitude of a sinusoidal signal. Right panels: pulsed 
fraction 3cr upper limits versus frequency in three different energy ranges (from top to bottom: 0.3— 2keV, 2— 4keV, 4-10 keV), calculated 
on the power spectrum binned at 0.00285 s (left panel) in the 0.005—175 Hz frequency range (see ^3. II for further details). 



3.3 Spectral analysis: continuum modeling and 
search for spectral lines 

The total spectrum, and the spectra relative to the 5 in- 
tervals shown in Fig.[3j were binned such to have at least 
60 counts per bin and not to over-sample the instrument 
spectral resolution by more than a factor of 3. We fit the 
total spectrum with an absorbed power-law (typical for this 
source; see e.g. Bosch- Ramon et al. 2007) using the XSPEC 
software, which gives a reduced xl = 1-33 for 173 degree of 
freedom (d.o.f.). The xi was n °t completely satisfactorily, 
we then studied the spectrum diving the total observation 
in 5 time intervals (see Fig. [3] and [6} to investigate whether 
this was due to a significant spectral variability within the 
observation. This was indeed the case. We fit together the 
spectra of the 5 time-intervals with an absorbed power-law 
function finding a, xt ~ 1-01 for 598 d.o.f, with the absorp- 
tion value tied to be the same for all 5 spectra and the slope 
and normalization of the power-laws free to vary. We find an 
absorption value of N H = (6.10 ± 0.08) x 10 21 cm~ 2 , and a 
power-law index which varied from 1.70 to 1.83 depending on 
the time-interval (see Fig. [3] all errors in the spectral values 
are reported at la confidence level). The source 0.5-10 keV 
absorbed flux varied from the first to the fifth time interval 
(from left to right in Fig.|3j): 8.5 ± 0.2, 6.6 ± 0.2, 7.2 ± 0.1, 
8.5±0.2 and 9.0±0.1 x 10" 12 ergs _1 cm" 2 (corresponding to 



unabsorbed fluxes of 1.2 ± 0.1, 0.9 ± 0.1, 1.0 ± 0.1, 1.2 ± 0.1 
and 1.3 ±0.1 x 10" 11 ergs~ 1 cm" 2 ). 

We also tried a fit using a pegpwr model instead of the canon- 
ical powerlaw model. The difference between the two is that 
in the former the normalization is calculated on an energy 
range given by hand (we assumed 0.3-10 keV), while in the 
latter case it is calculated at IkeV by default, and this can 
sometime cause a spurious F-normalization dependence. Us- 
ing a pegpwr we found T and flux values consistent within 
la with what previously reported for a simple power-law fit. 

We also tried to leave the absorption value free to vary 
among the 5 spectra. We found an absorption value vary- 
ing between 5.8(2)-6.2(2)xl0 21 cm~ 2 , with the same photon 
index or flux variation as the previous modeling. We then 
conclude that we do not detect any significant absorption 
variability within the sensitivity of our data. We do not 
find any evidence for absorption or emission features in the 
source averaged and time-resolved spectra. We derived up- 
per limits as a function of the line energy and width by 
adding Gaussian lines to the continuum model both in ab- 
sorption and in emission. For the total spectrum, the 3cr 
upper limit on the equivalent width of the presence of lines 
range between 50—550 eV depending on the energy (between 
0.3-10 keV) and assumed line width (cr E = 0.05 - 0.5 keV). 
As a simpler visual plot, we report in Fig.[S] (right panel) 
the la upper limit on a detectable flux variation with re- 
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spect to our continuum model as a function of the energy 
band, which corresponds to the detectability of spectral lines 
as broad as the energy-dependent instrument resolution on 
top of the continuum flux. 



4 DISCUSSION 

Thanks to this new Chandra observation (see we in- 
ferred the deepest upper limits to date on the presence 
of X-ray pulsations from the TeV Be binary LSI+61°303 
(PF < 10%; see 3370 and Fig.gJ, while the compact ob- 
ject was close to the apastron of its orbit around the Be 
companion (see Fig.[T|). Furthermore, we observed a strong 
flux variability during this observation (eventhough it was 
spanning a relatively small orbital phase in the 0.94-0.98 
range), in particular LSI+61°303 emitted two ks-timescale 
flares during our observation, with a rather hard spectrum 
with respect to the source persistent emission (see Fig.0 
and H3.3p . The accurate Chandra spatial resolution allowed 
us to study the X-ray spatial profile (although with the one- 
dimensional CC-mode capability; see H3.2p . which we found 
consistent with the instrumental PSF. 

4.1 X-ray pulsed fractions of rotational-powered 
and accreting pulsars, comparison with our 
current limits 

The detection of pulsations is the only unambiguous tracer 
for a secure determination of the pulsar nature of the com- 
pact objects hosted in LS 5039 and LS I +61°303. Deep 
searches for pulsations have been performed in the radio 
band at several frequencies, with the hope of detecting a 
fast spinning radio pulsar as in the case of the other TeV 
binary PSRB1259-63 (Johnston et al. 1999, 2005). How- 
ever, no radio pulsations have been detected so far from 
any of these two sources. This is not surprising; in fact the 
strong massive companion winds might have prevented ra- 
dio pulsations to be detected because of the strong free-free 
absorption all over the orbits and/or the large and highly 
variable Dispersion Measure (DM) induced by the wind. In 
particular, in that respect it is crucial to note that at peri- 
astron PSRB1259-63 does not show radio pulsations, and 
its periastron (given the large orbit, 3.4 year period) has 
about the same dimension of the major axis of the orbit 
of LSI+61°303 (26 days period) and it is way larger than 
LS5039's obit (4 days period), making very plausible the 
difficulty in detecting any radio pulsed emission from these 
closer binaries. 

On the other hand, searches for pulsations in the X- 
ray band have many more chances of success than in the 
radio band. In fact: the X-ray pulsar beam is usually larger 
than the radio one, and the strong companion wind does not 
influence much the X-ray pulsed emission if present. 

The deep PF upper limits we derived here (see [J3] and 
Fig. [4]) can be used to compare the X-ray emission of this 
TeV binary with that of isolated rotational-powered pulsars. 
X-ray pulsations from rotational powered pulsars have been 
detected in tens of objects, of which only fourteen are to 
date also TeV emitters (see e.g. Mattana et al. 2009), due to 
the presence of an energetic Pulsar Wind Nebula (PWN). 
All these TeV emitting pulsars known to date have high 
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Figure 5. Surface brightness of LS I +61°303. 1 pixel corresponds 
to 0.495". The inset represents the ratio between LSI+61°303 
and RX J0806+1527's surface brightnesses. 

rotational power (3 x 10 35 <£< 5 x 10 38 erg s" 1 ), fast 
pulsations (33 < P < 326 ms) , young characteristic ages 
(0.7 < r c < 90 kyr), rather high magnetic fields (10 12 < 
B < 5 x 10 13 Gauss), and X-ray pulsed fractions PF > 18%. 

The X-ray upper limits we derived here, namely 
PF ~10% (see Fig. 3} for spin periods as fast as ~ 5.6 ms 
(much faster than the fastest rotational powered pulsar 
known to be a TeV emitter), clearly indicates that if 
the compact object hosted in LSI+61°303 is indeed a 
rotational-powered pulsar, its X-ray emission cannot be 
dominated by synchrotron emission from the pulsar mag- 
netosphere alone, which might at most account for 10% 
of the emission, while the shock between the pulsar and the 
massive companion star winds is responsible for the remain- 
ing X-ray emission of the system. 

In fact in this pulsar scenario for LSI+61°303, radio, 
X-ray and gamma-ray emission are (at least in part) thought 
to come from the region where the pulsar wind interacts with 
the wind of the Be companion star and/or (particularly, the 
gamma-rays) from within the pulsar wind zone. This mostly 
depends on where the high-energy electrons generating the 
synchrotron and inverse Compton emission originate, the 
opacities therein, and the geometry of the orbit (Sierpowska- 
Bartosik & Torres 2008) . 

In this respect it should be mentioned that the TeV 
binary system hosting PSRB1259-63, does not show X-ray 
pulsations (Chernyakova et al. 2006) during its periastron 
passage, with a limiting pulsed fraction (derived as described 
in ^3. 1 P of ~ 15% (more details on that will be presented in 
the second paper of this series j3- Also for this source, the X- 
ray emission might in fact be dominated by the inter-wind 
and/or pulsar- wind zones rather than being due to emission 
from the pulsar magnetosphere. 

However, caution should be used here: in this line of rea- 
soning we are neglecting the possibility of a pulsar with a 

2 Note that Chernyakova et al. (2006) reports a limiting pulsed 
fraction of 2%, however these authors did not report the exact 
definition they used for pulsed fraction, the energy range, and 
background subtraction method they used. Possible differences 
in these assumptions might be the cause of such discrepancy. 



Deep Chandra observations of LS I +61° 303 7 







l| 1 INT-l: r=1.71(2) ^^^^^^ 
J| INT-2: T= 1.70(2) W 

INT-3: r=1.75(2) -f II . 

f" ' 

INT-5: r=1.83(2) 





Energy (keV) Energy (keV) 

Figure 6. Left panel: spectra of the 5 time intervals fitted with an absorbed power-law function (see Fig. \3\ and £|3,3I for details). Right 
panel: la sensitivity for the presence of spectral lines (defined as the la limit for the detection of flux variability as a function of energy). 



beam pointing away from our line of sight or spinning faster 
than ~ 5.6 ms, and also that strong pulsations might have a 
transient appearance in only a small fraction of the orbital 
motion, and in such small fraction being instead dominated 
by the pulsar magnetosphere. On the other hand, if we as- 
sume that a putative transient pulsed emission is diluted in 
part of the orbit by Compton scattering, from the rise times 
of the flares t ~ Rc t/c ~ 100 s (see Fig. [2] c is the velocity 
of light), and assuming an optical depth r >> 1, we can 
derive that the Compton scattering sphere radius should be 
Rc << 3x 10 12 cm, which is about half of semi-major axis of 
LS I +61°303's orbit. In this respect, our observations were 
performed in a part of the orbit (soon after the apastron 
passage; see Fig.[T]), where the detection of transient pulsa- 
tions should have been more favorable, and the Compton 
scattering highly reduced with respect to other parts of the 
orbit closer to the massive companion. 

Considering instead the accreting neutron star scenario 
(in the black hole case we obviously do not expect pulsa- 
tions at all hence we will not discuss this possibility which 
is of course open), low pulsed fractions are not unexpected. 
In particular, in many cases the pulsed X-ray emission com- 
ing from the accretion column might represent only a small 
fraction of the total X-ray emission which can also have a 
dominant un-pulsed component coming from the accretion 
disk and/or the neutron star surface. 

4.2 Flares, and flux/spectral variability on ks 
timescales 

A large variety of flares on several timescales (from seconds 
to days) have been observed to date from both neutron star 
and black hole binary systems. While in the accretion sce- 
nario flares are easily explainable as changes in the accretion 
rate from the companion star, or accretion disk, towards the 
compact object, this kind of X-ray variability has never been 
observed in isolated rotation powered-pulsarrjf). 

3 Note that magnetars do show X-ray flaring activity on several 
timescales, but in fact they are not rotational-powered while their 
emission (flaring and persistent) is due to the presence and insta- 



In the rotational-powered pulsar scenario, the ks- 
variability we observe here indicates that the interaction of 
the putative pulsar and the stellar wind change at distance 
scales which are smaller than the size of the binary orbit; 
put otherwise, that the mechanism by which the variability 
is generated has a shorter-than-orbital timescale. To inter- 
pret the short timescale variability (of the order of ks, and in 
the form of hard flares) seen in our Chandra data (also seen 
earlier in XMM-Newton observations by Sidoli et al. (2006) 
and Chernyakova et al. (2006)), it seems reasonable to in- 
voke dumpiness of the stellar wind (see ideas by Neronov & 
Chernyakova 2007, and Zdziarski, Neronov & Chernyakova 
2008). 

If one assumes that the power output in the pulsar wind 
is constant in ks timescales (which may not be necessarily 
true, due to inhomogeneities in the pulsar wind zone, and if 
so, ks variability could be related to it), the distance scale 
related to the fast variability has to be associated with the 
pulsar wind interactions with clumps of the Be stellar wind. 
In the latter case, the size, mass, and number density of the 
clumps imprint their signal on the X-ray emission. The dis- 
tance scale relevant for this fast variability can be estimated 
as: 




where v is typical velocity scale (e.g. the orbital velocity of 
the pulsar, or the speed of the stellar wind). 

The stellar wind of Be stars, like the one present in 
LS I +61°303, is typically assumed to have two components, 
one related with a polar contribution, and the other, with 
an equatorial disk (e.g. Waters et al. 1988). The polar wind 
is radiatively driven (e.g. Castor & Lamers 1979). The equa- 
torial wind is instead assumed to have an homogeneous ve- 
locity law V cq (r) = Vo q {r/R s ) m where m ~ 1.25, i? s is the 
stellar radius, and V e9 = 5kms _1 . These authors further 
assume that the terminal velocity of the equatorial wind 



bility of their ultra-strong magnetic fields (see Mereghetti 2008 
for a recent review). 



8 Rea et al. 



is V<^ ~ few 100 km s -1 , so that it is about a factor of at 
least 10 times smaller than the terminal velocity of the polar 
component, V^" ar . dumpiness in this wind can make the 
region of collision between the pulsar and stellar outflows to 
loose coherence, leading to the disappearance of a regular 
bow-shaped surface. If so happens, the escape timescale of 
the pulsar wind mixed with the stellar wind can be written 
as 

n/ inn ( ^ cm \ f W? cm M , 
tcsc ~ D/Vvrind ~ 100 I — I I — — — — I ks 

(slowed down from the case in which electrons move along 
the shock with the plasma velocity c/3, the drift velocity in 
the direction parallel to the shock surface), and where u w j n d 
can take typical or y^ lar values. 

Even with this timescale, clumps can cool electrons via 
inverse Compton. This statement can be entertained after 
comparing the cooling timescale of electrons of about 10 
MeV, which would be responsible via inverse Compton emis- 
sion of X-ray photons between 1 and 10 keV. Assuming that 
electrons diffuse in disordered magnetic field in the clump, 
one can estimate the Bohm-diffusion-regime escape time as: 

thus, clumps with such magnetic fields can retain electrons 
over significant timescales. The looses by synchrotron and 
inverse Compton are 

t synch ~ 3 x 10° (— j { - j s, 

, „ in4 ^ L * \ ( D \V l0MeV \ 
IC_ ^10 38 erg/sj VlO 12 cm) \ E ) S 

with L* being the luminosity of the star. For nominal val- 
ues of the parameters, then, an inverse Compton cooling is 
possible. 

Using the expression for the inverse Compton and syn- 
chrotron timescales one finds that if the magnetic field does 
not significantly rise from its pulsar wind nebula value of ~ 1 
G towards the center of the system, inverse Compton losses 
dominate at the inner region. However, synchrotron emission 
quickly dominate the farther out of the system electrons are, 
and can contribute to the overall X-ray emission detected, 
given the nebula an onion-like structure where losses dom- 
inance change. Depending on the location of the inhomo- 
geneity and value of magnetic field, relative contributions of 
inverse Compton and synchrotron X-rays can vary. 

Finally, if clumps can be described by variations in den- 
sity, and if the typical density of a clump in the wind is 
n c iump ~ 10 10 cm -3 and its size is 7? c i ump ~ 10 11 cm, then 
5Nh ~ n c i ump x -Rdump ~ 10 21 cm -2 . We did not find signif- 
icant ./Vh variability within our data, however the strong de- 
pendence between Nn, T and the power-law normalization, 
intrinsinc to the power-law modeling, might have hidden this 
variability in an increase of the hardness of the emission, as 
seen in our Fig. [3] and [6] (see also Smith et al. 2009). 



4.3 Limits on the presence of diffuse emission 

Studying the one-dimentional Chandra PSF, we did not find 
evidence for X-ray diffuse emission (see Fig. [5] and £I3.2[) . For 
a quantitative comparison of our limits with the hint for 
diffuse emission observed by Paredes et al. (2007; see their 
§2.3), we followed their analysis extracting the counts in 
an annular region of 5"-12.5"radii. Assuming an absorbed 
power law spectrum with Nh = 6.1 X 10 21 cm~ 2 and T = 2, 
we inferred a 3cr limiting unabsorbed flux for the presence 
of X-ray diffuse emission of ~ 1 x 10~ 14 erg s — 1 cm — 2 in the 
0.3-10 keV energy rang^f). Using the same assumptions we 
made, the flux of the diffuse emission proposed by Paredes 
et al. (2007) would be ~ 2 x 10" 14 erg s _1 cm" 2 , well above 
our detection limit. 

We caveat that if a putative extended emission (Paredes 
et al. 2007) shows a long term variability, our non detection 
can be due to a different observing epoch with respect to the 
original Chandra observation which hinted at the presence 
of diffuse emission. 



5 SUMMARY 

With the Chandra observation we report here, we in- 
ferred the deepest upper limits to date on the presence 
of X-ray pulsations from the TeV Be binary LSI+61°303 
(PF £j 10%; with a strong energy and frequency depen- 
dence, see Fig. [4]), while the compact object was close to the 
apastron of its orbit around the Be companion (see Fig.Q]). 
These PF limits are deeper than the pulsed fraction mea- 
sured from any isolated rotational-powered pulsar, in par- 
ticular from those emitting in the TeV range, while they are 
well in line with what is observed from accreting neutron star 
binaries (although we do not see clear signs of accretion, as 
e.g. spectral lines). Hence, if hosting a rotational-powered 
pulsar, LS I +61°303's X-ray emission cannot be driven by 
the same processes as for a similar isolated case (with the 
caveats discussed in ^4. II though), while indeed a major con- 
tribution to its X-ray emission should come from the shock 
between the pulsar and Be star winds. This is probably also 
the case for PSRB1259-63. The ks-timescale flares we saw 
in this Chandra observation of LSI +61° 303 show a harder 
spectrum with respect to the source emission, and can be 
explained in an accretion scenario by variability in the ac- 
cretion rate, while in the rotational-powered pulsar one, as 
the interaction between the pulsar wind and clumps in the 
Be wind. 

This research has made use of data from the Chandra 
X-ray Observatory and software provided by the Chandra X- 
ray Center. NR is supported by a Ramon y Cajal Research 
Fellowship to CSIC, and thanks G.L. Israel for useful discus- 
sion, advices concerning the timing analysis, and for allowing 
the use of his DPS software. We thank the anonymous referee 
for his/her very useful suggestions, and V. Bosch-Ramon 



4 Note that to infer the diffuse flux, Paredes et al. (2007) assumed 
the spectrum of their Chandra observation, T = 1.25, which was 
highly affected by pile-up. We then prefer to assume T = 2, the 
typical value for a PWNe. 
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and M. Ribo for comments. This work has been supported 
by grants AYA2009-07391 and SGR2009-811. 



REFERENCES 

Abdo A. et al. 2009a, ApJ Letters 701, 123 

Abdo A. et al. 2009b, ApJ Letters 706, 56 

Aharonian F., et al. 2005a, A&A, 442, 1 

Aharonian F., et al. 2005b, Science 309, 746 

Aharonian F., et al. 2006, A&A 460, 743 

Aharonian, F., et al. 2007, A&A 469, LI 

Albert J. et al. 2006, Science 312, 177 

Albert J. et al. 2007, ApJ 665, L51 

Albert J. et al. 2008, ApJ 684, 1351 

Albert J. et al. 2009, ApJ 693, 303 

Aragona, C, et al. 2009, ApJ 698, 514 

Barthelmy, S., et al. 2008, GCN 8215 

Bosch-Ramon, V., et al. 2005, A&A 628, 388 

Bosch-Ramon, V., et al. 2007, A&A 473, 545 

Casares, J., et al. 2005, MNRAS, 364, 899 

Castor J. I., & Lamers H. J. G. L. M., 1979, ApJ, 39, 481 

Chcrnyakova, M., Lutovinov, A., Rodriguez, J., Revnivtsev M., 

2005, MNRAS, 364, 455 
Chernyakova, M. et al. 2006, MNRAS, 372, 1585 
Chcrnyakova, M., Neronov, A., Aharonian, F., Uchiyama, Y., 

Takahashi, T. 2009, MNRAS, 397, 2123 
Clark, J., et al. 2002, A&A, 392, 909 

Cottam, J. Sako, M., Kahn, S. M., Paerels, F., Liedahl, D. A., 

2001, ApJ, 557, L104 
de Pasquale, M. et al. 2008, GCN 8209 
Dhawan, et al. 2006, Micro, work. p. 52 
Dubus G. 2006, A&A 451, 9 
Dubus, G. & Giebels 2008, Atel#1715 

Esposito, P., Caraveo, P. A., Pellizzoni, A., de Luca, A., Gehrels, 

N., Marelli, M. A. 2007, A&A 474, 575 
Frank, King & Raine 1992, Science 258,1015 
Gregory, P.C. 2002, ApJ, 525, 427 
Grundstrom, E.D. et al. 2007, ApJ, 656, 437 
Harrison, F. A., et al. 2000, ApJ 528, 454 

Hickox R. C, Narayan R., Kallman T. R, 2004, ApJ 614, 881 
Hinton, J., et al. 2009, ApJ 690, L101 
Israel, G.L. & Stella, L. 1996, ApJ, 468, 369 
Israel, G.L., et al. 2003, ApJ, 598, 492 

Johnston, S., Manchester, R.N., Lyne, A., Bailes, M., Kaspi, V. 
M., Qiao, G., DAmico N., 1992, ApJ ,387, L37 

Johnston, S. Manchester, R.N., McConnell, D., Campbell- Wilson 
D., 1999, MNRAS, 302, 277 

Johnston, S., Ball, L.,Wang, N., Manchester, R. N., 2005, MN- 
RAS, 358, 1069 

Kishishita, T., Tanaka, T., Uchiyama, Y., Takahashi, T. 2009, 

ApJ, 697, LI 
Leahy, D. A., 1987, A&A, 180, 275 

Maier G. et al. 2009, Proc. of the 31st ICRC, Lodz 2009, 0512 
Massi M., et al. 2001, A&A 376, 217, 
Massi M., et al. 2004, A&A 414, LI 
Mattana F., et al. 2009, ApJ 694, 12 

Munoz-Arjonilla et al. 2008, Atel #1740 

Neronov A, & Chernyakova M. 2007, |astro-ph/0701144l 

Paredes, J. M., et al. 2007, ApJ 664, L39 

Ray, P. et al. 2008, Atel#1730 

Rea, N. & Torres, D.F. 2008, Atel#1731 

Ribo, M., Paredes, J. M., Moldon, J. Marti, J., & Massi, M. 2008, 

A&A, 481, 17 
Sidoli, L. et al. 2006, A&A 459, 901 

Sierpowska-Bartosik A. & Torres D. F. 2008, Astroparticlc 

Physics 30, 239 
Sierpowska-Bartosik A. & Torres D. F. 2009, ApJ 693, 1462 



Smith, A. W., Kaaret, P., Holder, J., Falcone, A., Maier, G., Pan- 
del, D., Stroh, M., 2009, ApJ 693, 1621 

Torres, D. F., Sierpowska-Bartosik, A. & Rea, N. 2010, Reports 
on Progress in Physics, submitted 

van der Klis, M. 1989, Timing Neutron Stars, Eds. H. Ogelman 
& E.P.J, van den Heuvel (Dordrecht: Kluwer), 27 

van der Meer, A., Kaper, L., di Salvo, T., Mendez, M., van der 
Klis, M., Barr, P., Trams, N. R., A&A, 432, 999 

Vaughan, B. A., et al. 1994, ApJ, 435, 362 

Waters L. B. F. M., et al., 1988, A&A, 198, 200 

Zdziarski A. A., Neronov A, & Chernyakova M. 2008, 
larXiv:0802.1174l 



